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Research Objectives 

The objectives for this SOLVE project were 3-fold. First, we sought to calculate a complete set 
of photolysis rate coefficients (j-v alues) for the campaign along the ER-2 and DC-8 flight tracks. 
En route to this goal, it would be necessary to develop a comprehensive set of input geophysical 
conditions (e.g., ozone profiles), derived from various climatological, aircraft, and remotely sensed 
datasets, in order to model the radiative transfer of the atmosphere accurately. These y- values 
would then need validation by comparison with flux-derived y-value measurements. The second 
objective was to analyze chemistry along back trajectories using the NASA/Goddard chemistry 
trajectory model initialized with measurements of trace atmospheric constituents. This modeling 
effort would provide insight into the completeness of current measurements and the chemistry of 
Arctic wintertime ozone loss. Finally, we sought to coordinate stellar occultation measurements 
of ozone (and thus ozone loss) during SOLVE using the MSX/UVISI satellite instrument. Such 
measurements would determine ozone loss during the Arctic polar night and represent the first 
significant science application of space-based stellar occultation in the Earth’s atmosphere. 


Summary of Progress and Results 

Photolysis Rate Coefficients during SOLVE 

The model calculation of 7-values plays a vital role in the photochemical analysis of the chem- 
istry of ozone production and loss in the SOLVE campaign. Building upon our 7-value modeling 
success in POLARIS [Swartz et al., 1999], the sophistication of the input of geophysical data to 
the APL radiative transfer model was increased in order to cope with the high-solar zenith angle 
(SZA) conditions found in polar winter. A comprehensive ozone/pressure/temperature climatology 
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was adopted as a baseline description of the atmosphere. Earth Probe Total Ozone Mapping Spec- 
trometer (TOMS) satellite and Composition and Photodissociative Flux Measurement (CPFM) 
spectroradiometer [McElroy, 1995] data were utilized, as in POLARIS, along with the addition of 
remote sensing data from the Polar Ozone and Aerosol Measurement (POAM) III satellite instru- 
ment and in situ ozone data. The model inputs have been compared and the impact on 7-value 
calculations examined along the ER-2 and DC-8 flight tracks. 

The DC-8 /-values were also compared with the Scanning Actinic Flux Spectroradiometer 
(S AFS) [Shetter and Muller , 1 999], which measured the An actinic flux on board the DC-8, for the 
first time. The overall model-measurement agreement for the DC-8 was excellent, within ~5% 
for 7 no 2 ( n0 2 -> NO + 0( 3 P)) and ~ 10% for j 0j (0 3 -y 0 2 + 0( J D)). This extensive compar- 
ison also revealed a curious bias in Jno 2 that is consistent with the attentuation of light by polar 
stratospheric clouds (PSCs). In Figure 1, model calculations utilizing TOMS (APLtoms) and cli- 
matological (APL clim ) data are compared with SAFS measurements. In the 20000123, 20000125, 
and 20000308 flights, where PSCs were particularly prevelent, model— measurement agreement de- 
viates significantly from the average. During these high-PSC flights, the measured 7- values were 
lower than the model calculations. The effect is most noticable at high SZAs (Figure lb), when 
Sunlight would have passed nearly horizontally through the putative cloud layers. 


NASA/Goddard Chemistry Trajectory Modeling of SOLVE in situ Measurements 

The NASA/Goddard trajectory chemistry model, initialized with ER-2 data, has been used to eval- 
uate heterogeneous processing and radical photochemistry within the winter 1999—2000 Arctic 
polar vortex. From this preliminary look, we have found that the model reproduces midlatitude 
and polar extravortex chemistry with a level of model-measurement agreement comparable to that 
found during POLARIS [Pierson et al. , 2000]. Within the vortex, however, there are inconsisten- 
cies that warrant further study. 

HQ* partitioning. Similar to POLARIS, the model under-predicts extravortex OH and H0 2 in- 
dividually by about 20% and OH/HO2 by approximately 5% on average. The model over-predicts 
OH/H0 2 by more than 50% within the vortex, using JPL97 cross sections/quantum yields and ki- 
netics, with the addition of more recent OH+NO, data. The use of JPL00 7-data leads to better 
agreement in H0 2 outside the vortex (see Figure 2, for example), with OH/H0 2 still matching 
well. Using JPL00 k - data and 7-data increases H0 2 outside the vortex (now over-predicting H0 2 ) 
and worsens the OH/H0 2 agreement. OH/H0 2 (when the Sun is up) is dominated by H0 2 +C10, 
OH+O3, and OH+CIO (HOCl+hv is often marginally rate-limiting). OH+CIO must be included 
in the steady-state HO x equation inside the vortex. Steady-state calculations based on model con- 
centrations agree well with modeled values (good approximation of important HO x chemistry). 
Extravortex HO x steady-state calculations using measured NO, 0 3 , CO, and CIO are lower than 
measured OH/H0 2 and about 50% high within vortex. This is slightly improved by subsituting the 
H0 2 +C10 term withy Hocl /H0 2 , dependent on the choice of measured or modeled [H0 2 ], HO x 
within the vortex is improved further by using Cl/ (below). 


HC1 measurements. During many of the SOLVE vortex survey flights, the model over-predicts 
HC1 destruction (even with heterogeneous processing disabled). This disagreement disappears in 
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Figure 1 . Comparison of (a) APL X0MS and SAFS (measured; SZA < 85°) and (b) APL clim and 
SAFS (measured; SZA > 85°) y N Q 2 along the DC-8 flight tracks during SOLVE. Each science 
flight is identified by a unique symbol, as keyed at right. 
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SOLVE 2000031 1 Trajectory Chemistry 



Figure 2. Inside-outside polar vortex comparison of in situ measured trace gases with modeled 
chemistry along trajectories for the March 11, 2000 SOLVE ER-2 flight. Measured quantities 
(solid lines) are shown, along with calculations made with the Goddard trajectory chemistry model 
along 10-day trajectories (+) using DAO meteorological fields. The chemical model was initialized 
with measured NO*, NO y , CIO, C10N0 2 , HC1, Cly, Br y , HO*, H 2 0, in situ 0 3 (fixed), overhead 
0 3 column, aerosols, and surface albedo. The shaded area indicates the extra-vortex segment of 
the flight. The accuracy of modeling of the HO* budget outside the vortex is comparable (or better) 
to that found during POLARIS. Significant discrepancies are found within the vortex, however. 
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March as the vortex warmed and weakened. Temperatures must be increased uniformly by 2-5 K 
to keep prevent HC1 from being destroyed. HC1 measurements are not correlated with solar zenith 
angle, latitude, N 2 0, temperature, or altitude (but there is greater HCI variability at around 20 km). 


Cly partitioning, Cly is in good agreement with the N 2 0 correlations of R. J. Salawitch. There 
is too much Cl activation (CIO, C10N0 2 ) within the vortex. Using JPL00 data provides no CIO* 
improvement. Use of Cl/=measured HCl+C10N0 2 +CIO+2C100Cl/non-(HOCl+2Cl 2 ) fraction 
of modeled C\ y reduces Cl activation. Cly*, however, does not equal Cly, even outside the vortex, 

so there is likely a problem with the measured Cly budget. 

Photochemical Ozone Loss Determined by MSX/UVISI Stellar Occultation 

The combined SOLVE campaign sought, in part, to quantify ozone loss within the Arctic po- 
lar vortex using a variety of aircraft-, balloon-, ground-, and space-based instrument platforms. 
The Midcourse Space Experiment (MSX)/Ultraviolet and Visible Imagers and Spectrographic Im- 
agers (UVISI) suite of instruments performed 3 1 stellar occultation observations from January 23 
through March 4, 2000, in and near the Arctic polar vortex. Using a newly developed, combined 
extinctive-refractive algorithm [Yee et al., 2001], ozone mixing ratio profiles, along with profiles of 
total density, pressure, and temperature, were retrieved. Retrieved temperature and ozone profiles 
were shown to agree well with other measurements. Diabatic trajectory calculations were used to 
remove the effects of subsidence within the vortex, allowing photochemical ozone loss to be in- 
ferred from these occultation measurements. A maximum ozone loss of about 1 ppmv was found 
at 400-500 K (—16-21 km), as shown in Figure 3, during the 41-day period for which occultation 
data are available, in agreement with several other ozone loss analyses of the campaign [Newman 
et al., 2001]. This corresponds to an average daily loss rate of -0.024 ppmv/day. If we assume 
an average loss rate of 0.05 ppmv/day for the second week in March, a total loss of 1.4 ppmv for 
1/20-3/12/2000 is inferred, corresponding to the comparison reported by Newman et al. [2001], 
The combined extinctive-refractive stellar occultation technique has now been demonstrated to 
accurately measure stratospheric ozone loss during polar winter [Swartz et al., 2001]. 



6 


5. A. Lloyd and W. H. Swartz 


Average Rate of Ozone Change (ppmv/day) 
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Figure 3. Profiles of the average daily rates of ozone change and total ozone destruction, from 
January 23 to March 4, 2000. The results from three different but related methods, utilizing average 
vortex descent rates and two analyses using individual diabatic trajectory calculations, are shown. 
The data have been averaged over 25-K layers, and the lined region contains the average l-cr 
uncertainties of the linear fits of the ozone measurements (and thus the rates of ozone change) on 
the individual diabatic surfaces within each layer. 
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